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ABSTRACT 

Transiting cxoplanetary systems are surpassingly important among the planetary sys- 
tems since they provide the widest spectrum of information for both the planet and 
the host star. If a transiting planet is on an eccentric orbit, the duration of transits 
Td is sensitive to the orientation of the orbital ellipse relative to the line of sight. 
The precession of the orbit results in a systematic variation in both the duration of 
individual transit events and the observed period between successive transits, P bs- 
The periastron of the ellipse slowly precesses due to general relativity and possibly 
the presence of other planets in the system. This secular precession can be detected 
through the long-term change in P bs (transit timing variations, TTV) or in Td (tran- 
sit duration variations, TDV). We estimate the corresponding precession measurement 
precision for repeated future observations of the known eccentric transiting exoplane- 
tary systems (XO-3b, HD 147506b, GJ 436b and HD 17156b) using existing or planned 
space-borne instruments. The TDV measurement improves the precession detection 
sensitivity by orders of magnitude over the TTV measurement. We find that TDV 
measurements over a ~ 4 year period can typically detect the precession rate to a 
precision well exceeding the level predicted by general relativity. 

Key words: binaries: eclipsing - planetary systems - relativity - methods: observa- 
tional - techniques: photometric 



1 INTRODUCTION 

Since the discovery of the first transiting ext rasolar planet 
l|Charbonneau et al.l l2000l : iBrown et al.l 1200 ll ). the number 
of such systems has increased to more than 3Qj. These tran- 
siting extrasolar planets (TEPs) provide unique information 
on the properties of the system. Based on the geometry pro- 
vided by the transit light curve(s), the inclination, the phys- 
ical radius and mass, therefore the density and the surface 
gravity can be derived, in addition to the mass of the planet. 
Moreover, the time between successive transits can be mea- 
sured with an exceedingly high accuracy (~ 10 -6 - 10 -7 , 
relative to the period). The detection of long-term tran- 
sit timing variations can be used to l earn more beyond the 
properties of the par ent-star system (|Miralda-Escudeil2002l : 
ISteffen fc AgoT2 007). They can be indicative of the presence 
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of other planetary companions (see e.g. iHolman fc Murravl 
120051 : lAgol et all 120051 : iMiller-Ricci et all l2008h . co-orbital 



comp anions (Trojans, s ee iFord fc Holmanll2007l ). or satel- 
lites (|Simon et al.l 120071 ) in the system, could provide in- 
formation on the oblateness of the host star, or can be 
used to detect the additional prograd e periastron precession 
predicted by general re lativity (GR) l|Miralda-Escudell2002l : 
iHevl fc Gladmanll2007l ). Secular variations in the semimajor 
axis (and therefore in the transit timing) are also predicted 
on the time scale of stellar life due to the anisotropic ligh t 
redistributio n (a.k . a. Yar kovski- effect, see iFabrvckvl 120081 ) ■ 
Furthermore. Iloriol (|2006h has shown that TEP observations 
can in principle also test the gravitoelectric correction of 
GR by measuring the radial velocity amplitude and tran- 
siting periodicity simultaneously, in order to verify that the 
third Kepler's law requires a semimajor axis dependent cor- 
rection. 



In a pioneer study, iMiralda-Escudd (|2002l ) derived the 
modification of the observed time period between succes- 
sive transits P bs, called transit timing variations (TTVs), 
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caused by the standard periastron precession due to GR (e.g. 
M isner. Thorne. fc Wheeleilll973l ) and the perturbat i ons o f 
other planets if present. Recently, iHevl fc Gladmanl (2007) 
have extended these studies and estimated the precision of 
precession rate measurements for long-term mock obser- 
vations of eccentric transiting extrasolar planets (ETEPs). 
Both studies restricted to small eccentricities. At that time, 
the existence of close eccentric planets was known only 
through radial velocity measurements, and no ETEPs had 
been observed. Since their publication, four transiting ex- 
trasolar planets have been discovered w ith significant eccen- 
tricity: X O-3b llJohns-Krull et ai]l2007t ). HP 147506b (a.k.a 
HAT-P-2. iBakos et al.l l 2007l ,, GJ 436 b dGillon et all 120071 ; 
iButler et al.ll2004l ). and HD 17156b l|Fischer et al l l2007f ) . 
Therefore it is now possible, for the first time, to make 
specific predictions for future, long-term measurements of 
periastron precession effects for real exoplanetary systems. 

In this paper we determine the precision by which re- 
peated long-term future ETEP observations will be able to 
detect the periastron precession rate for existing systems. 
In addition to TT Vs, i.e. the slow modu l ation of P b s con- 
sider ed previously l)Miralda-Escudel 120021 ; IHevl fc Gladmanl 
120071 ) , the periastron precession also changes the time dura- 
tions Td of individual transits. We examine whether these 
transit duration variations (TDVs) can be used to improve 
the sensitivity of periastron precession measurements. We 
estimate the precession rate measurement precision for long 
term repeated observations of P bs and Td for the known 
ETEPs. Since several of the observed ETEPs have large ec- 
centricities, we derive expressions for both TTVs and TDVs 
which are applicable for arbitrary eccentricities. We estimate 
whether future observations of currently known ETEPs will 
be able to reach the sensitivity necessary to test the pre- 
diction of GR, using existing or planned space-borne instru- 
me nts. We refer the reade r to a recent independent study 
by I Jordan fc Bakoj (|200Sf ). of precession rates in eccentric 
transiting extrasolar planets. 

The next section of this paper introduces the geometri- 
cal description which is the basis of our calculations, and de- 
rives the expected transit timings and durations for planets 
orbiting a star with an arbitrarily large eccentricity. In § O 
we utilize our results for the confirmed four ETEP systems, 
and give predictions for future observations of periastron 
precession with space-borne observations. Our conclusions 
are discussed in § 2] 



2 TRANSIT TIMINGS AND DURATIONS FOR 
ECCENTRIC ORBITS 

The reference frame used for the description of exoplane- 
tary systems as well as for binary/multiple stellar systems 
is fixed to the sky: the plane of the sky is defined by the 
(X+, Y+) while Y+ points towards to north. For planetary 
transit observations, the line-of-sight lies close to the orbital 
plane, i.e. perpendicular to the plane of the sky. The orbit is 
given by Cartesian coordinates (£, rf), where £+ is parallel to 
V— and 77+ oriented toward the observer (see also Fig. [T|. 
The Lagrange vector or eccentricity vector is given in these 
coordinates as (fc, /i) = (ecostj, esintj). Let us define the 
angle ipo as the angle relative to £+ in the orbital plane at 




(observer, Earth) 

Figure 1. The geometry of the orbit of the transiting planet. The 
plane of the sky is defined by the X+ (west) and Y+ (north) axes 
while the Z+ axis points away from the Earth. The orbital plane 
is defined by the axes £+ and rj+ , where £+ is anti-parallel with 
X+ axis, r?+ is in the plane of (Y+, Z+) and the angle between 
Y+ and rj+ is the inclination, nearly 90°. The major axis of the 
orbit is marked by the dotted line. 



the instancfl of t fle transit. From the definitions of (£,77), 
observing from Earth is equivalent to setting ipo — tt/2. 

Now let us denote the mean longitude of the planet at 
the instance of the transit by A. For a circular orbit, A = 
ipo. From its standard definition in celestial mechanics, it is 
straightforward to derive the mean longitude for arbitrary 
eccentricities (see Appendix [XJ. The result is 

A = X(ipo, fc, h) = \(tfio, e cosu;, e sinoj) = 

/ he± . fcej 

= arg I fc + cos ip + -j—j , h + sin ipo - 75— 

ejfl-l) 
1 + eii ' 



(1) 



where en = fccosipo + h sin ipo, e± = fcsinyo — hcos<fio are 
the components of the eccentricity vector relative to the line- 
of-sight, I = 1 — vl — e 2 is the oblateness of the orbit, and 
arg(x,y) = arctan(y/s) if x ^ and 7r + arctan(y/x) other- 
wise. Plugging in the observed values of e and to, equation |T} 
provides a simple way of calculating the mean longitude of 
the orbit for an arbitrary transit observation. Note that this 
formalism omits the direct usage of the mean anomaly, true 
anomaly and eccentric anomaly which have no meanings for 
e — > 0. All of our derived formulae are based on the well- 
behaved parameters A, fc, and h, and therefore are valid for 
arbitrary eccentricities. 

In the following subsections, we derive the expressions 
describing TTVs and TDVs, discuss the corresponding ob- 
servational implications and estimate the precision of peri- 
astron precession observations. 



2.1 Modulation of the transit period 

The period between successive transits P bs is modified by a 
possible slow precession of the orbital elements. These mod- 
ulations are referred to as transit timing variations (TTVs) . 



2 i.e. at the center of a transit 
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Table 1. Basic data of the four known eccentric transiting exoplanetary systems: the mass (M*) of the parent star, period (P, in days), 
eccentricity (e) and the argument of pericenter (ui), the half-duration of a transit event (H, in days) and the impact parameter (b). In the 
last two columns, we provide the calculated values of the secular period caused by the GR periastron precession, and the minimum mass 
- semimajor axis ratio for a hypothetical exterior perturber at a? > 3a that would lead to the same magnitude of periastron precession. 



System 


M*/M 


P(d) 


e 


u> (degrees) 


H (d) 


b 


Psec (years) 


m 2 /M 9 
(a 2 /3a)3 


HD 147506b 


1.298 ±0.07 


5.63341 


0.520 ±0.010 


179.3 ±3.6 


0.083 





19790 ± 740 


12.2 


XO-3b 


1.41 ± 0.08 


3.19154 


0.260 ± 0.017 


344.6 ± 6.6 


0.050 


0.8 


9280 ± 410 


15.9 


GJ 436b 


0.452 ±0.013 


2.64385 


0.150 ±0.012 


351.0 ± 1.2 


0.065 


0.85 ±0.02 


15180 ± 400 


2.6 


HD 17156b 


1.2±0.1 


21.21725 


0.6717 ± 0.0027 


121.23 ± 0.40 


0.098 


0.50 ±0.12 


143000 ± 7900 


5.9 



We derive the modulation of P bs due to periastron preces- 
sion in two steps. First we demonstrate that the change in 
the period between successive transits is simply related to 
the change in the mean longitude AA. Then relating the 
mean longitude shift to the periastron precession rate we 
derive the expected TTV rate. 

Let Po be the orbital period, n = 2it/Po be the mean 
angular velocity. The mean longitude of the planet increases 
steadily in time, A = nt + Xq. The variation in the mean lon- 
gitude at the transit center would result in a variation in the 
transit cadence. During a transit at time ti, the mean longi- 
tude is Atr(ii) = nti + Ao, and after an observed revolution, 
at the instance ti = ti + Pobs it is At r (fa) + 2tt — nti + A u . 
Therefore the observed period between transits is 



2n + AA 

Pobs = = Po 



AA 
2^T 



(2) 



where AA = A tr (i2) — A tr (ii). 

In the following we assume that the shift in the mean 
longitude is caused by the perihelion shift Auj = Pquj per 
period, (i.e. we assume a constant eccentricity), then from 
the chain rule 



OX . OX Ok OX Oh . 

Olu Ok Ouj Oh Ouj 



(3) 



Here 0k /Ouj — — esinu; and Oh/Ouj = ecosw, from definition 
(see above), and the partial derivatives OX /Ok and OX /Oh 
can be found from equation ((TJ and are given explicitly in 
the Appendix (|A8HA9|) . At transit, we get 



OX = e 2 
0lo ~ 1 + VI - e 2 



-e 2 + (2 + esintj)esina; 



(4) 



(1 + esinu;) 2 

Combining equations ©-Q and defining the secular period 
of the periastron precession as P sec = (2tv)/u, we get 



P _ P ,Po 2 9X OX Pi 

Jobs — P} + 7T~ — P> + T{ „ • 

2ty Ouj Ouj Psec 



(5) 



Since OX/ Ouj itself is not constant due to periastron pre- 
cession, the observed period between transits slowly changes 
Differentiating equation © with respect to time, we get 

<9 2 A P 2 _ 47r(l-e 2 ) 3/2 ecosw P 2 



Pobs = 2-7T- 



(6) 



0uj2Pg c (l±esina;)3 P 2 e ' 

since the partial derivative of equation @ with respect to 
uj is 

2 X _ 2(1 - e 2 ) 3/2 ecosu; 
Ouj 2 (1 + esino;) 3 



(7) 



Note that this equation clearly shows that the small ec- 
centricity approximation 2 X/0uj 2 ~ ecosw (e.g. used by 



iMiralda-Escudel |2002| ; iHevl fc Gladmanl I2007T I is very im- 
precise for moderate to large eccentricites. Depending on 
the actual value of uj, equation (0 can result even 6 — 8 
times smaller or larger values for 2 X/Quj 2 as its first order 
approximation for eccentricities 0.5 — 0.7. 

We have to mention here that the observed period and 
therefore the individual transit timings are also affected by 
the light time effect (LTE). Since the precession of an ellip- 
tical orbit causes the distance between the host star and the 
planet at the transit instances to vary, the light time delay 
will change for transit event to transit event. The magni- 
tude of this effect can be derived as follows. The distance 
between the star and the planet at transit time (i.e. when 
X — (po = 7r/2) is 



a(l-e 2 ) 

1 + ecosu 1 + e cos(</?o — uj) 



a(l-e 2 ) 



a(l-e 2 ) 



(8) 



This difference in the distance implies an additional — r(l — 
fi)/c time shift, respective to the barycentric reference frame 
of the star-planet system. Here, fj, = M P /(M P + M*) <§C 1 
is the mass parameter, a is the semimajor axis of the orbit 
and c is the speed of light. Therefore, the correction in the 
observed period is 



3+LTE 



= P„> 



= P„, 



a n 2 ,0(l + esmuj) 1 . 

-(1 - e ) Pquj - 

c Ouj 



2,^ 



e )ecosw Pr 







c (1 + esinix)) 2 P sc 



0) 



where we neglected the barycentric correction. Thus, the 
variation in this period (corrected for LTE) due to the vari- 
ation in uj is 



fots™ = Pobs-4^ 2 (l-e 



2 e(e + ecos 2 uj + sinuj) a Po 



(1 + esinoj) 3 



(10) 



Comparing equation © and (|10p . and assuming that the 
motion of the planet is non-relativistic, a/c Po, we find 
that |P ts TB - ^obs| < I Pobs |- We conclude that the period 
variation due to the LTE is negligible compared to the period 
variation caused by the changing geometry. 

In summary, equations (0)-(j6]), along with equation Q, 
give the modulation of the actual observable period between 
transit events. These results are valid for arbitrary eccentric- 
ities and are independent of the physical mechanism causing 
the secular precession of the periastron. We calculate the 
secu lar precessi on period caused by general relativity (see 
e.g.. lWald1ll984l ) using 



p — 

i sec — 



c 2 (l-e 2 ) p5/3 
3(27rGM*) 2 / 3 ' 



(11) 
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where is the mass of the parent star, G is Newton's 
gravitational constant. This secular period is of order 10 4 - 
10 5 years for the known ETEP systems (see Table [1] and 
Section [3] for more details). We note that if other planets 
are also present in these systems, they might also cause ad- 
ditional periastr on precession of a larger magnitude. The 
Yarkovs ki-effect (|Fa brvckv 20081) and the tidal circulariza- 
tion (see lJohns-Krull et al.l2007l , and the references therein) 
lead to negligible modifications for our purposes, as these ef- 
fects are relevant on timescales exceeding 0.1 Gyr for these 
systems. In the following, we compare the precession mea- 
surement sensitivities with the general relativistic rate P aec 
given by equation ([lip. 



2.2 Modulation of the transit duration 

Here we investigate how the periastron precession affects the 
duration of a transit. Let us denote the half duration of the 
transit by H = Tb/2, which we define as half the time be- 
tween the instances when the center of the planetary disk 
intersects the limb of the star, i.e. between the center of the 
ingress and egress. Note that this is not the time between 
the first and last contact. This is important because the 
instances of the center of the ingress and egress can be mea- 
sured more accurately than their beginning or end. Since we 
are interested in estimating the variations of the duration of 
the transit to leading order, we perform a first-order calcu- 
lation, i.e. assuming a constant apparent tangential velocity 
for the planet and neglecting the changes in the impact pa- 
rameter due to the elliptical orbit and/or the curvature of 
the projection of the orbit due to the inclination. From Ke- 
pler's Second Law and equation JH}, one can calculate the 
tangential distance Aa: traveled by the planet during time 
H, 



1 + esinw 
Ax = Utan-H = an — H 



(12) 



(see Appendix iBl for the derivation of wtan). This can be 
related to the impact parameter b and the radius of the star 
Ri, for the geometry of the transit as 

^ = yr^. (is) 

Thus, to leading order, 



H = 



P I P* VT~~e5 



2tt I a 1 + e sin u 



y/l-b 2 + 



(14) 



Note that equation (|14[) depends on lu through the (1 + 
esinu;) -1 term and also implicitly through b, 



b = 



P* 1 + e sin lu 



(15) 



The variation in H caused by the variation in the periastron 
can be found from equation (|14[l and l)15[). 



dH _ _ecoso^^l_-2P 
8lu 1 + e sin to 1 — b 2 



(16) 



Note that equation reflects the qualitative expections im- 
plied by Kepler's Second Law. Namely, if an eccentric orbit 
advances, the distance between the planet and the star will 
change. If this distance decreases, the impact parameter will 
also decrease (resulting a longer transit duration) but due 



to Kepler's Second Law, the apparent tangential velocity of 
the transiting object will increase (resulting a shorter tran- 
sit duration). Therefore at a certain value of the inclination 
and/or the impact parameter, the two effects cancel each 
other yielding no TDV effect. Equation[l6]clearly shows that 
it occurs when the impact parameter is b — l/y/2 « 0.707. 
The long-term variation in the duration of the transit is then 



• _ dH _ VI - e 2 ecosu 1 - 2b 2 P* P 
~ ~d^ U ~ (1 + esincj) 2 ^1 - b 2 ~a~R^ 



(17) 



Comparing equation ([6} and equation (|17[) the TDV effect 
relates to the TTV effect as 



H _ 1 + esinw 1 - 2b 2 R* 

Rohs 



t>7T 



VT - b 2 Rsch 



(18) 



where Rsa = 2GM/c 2 is the Schwarzschild radius of the 
star. As an example, note that P*/Ps c h = 2.5 x 10 for 
the Sun. Therefore, equation (|18[l shows that the TDV ef- 
fect is always much larger than the TTV effect. In partic- 
ular, the ratio is larger for increasing b. In the limit 6 — *• 1 
equation (|17[) breaks down because the periastron precession 
shifts the orbit out of the transiting region. 



2.3 Observational implications 

Now, using the results for the TDV and TTV effects, 
equation © and (|17[) . we can estimate how these timing 
and transit duration variations might be observed on long 
timescales. In the following we discuss these observational 
implications. 



2.3.1 Transit Timing Variations 

Equation (J6| shows that the observed period between suc- 
cessive transits increases or decreases at a practically con- 
stant rate during the observations, P bs- The time of the 
mth transit from an arbitrary epoch To can be found from 
adding up the contributions of the observed m number of 
periods 



T + P ohs m + Dm^ 



(19) 



where P b s « Po denotes the time of the first observed orbit, 
D is the transit timing variation factor, 



,-, PoPobs „ „ / . / Po 

D = — - — = 27rG T v(e,u)P — 

^ \ Jscc 



where we have introduced the geometrical factor 



GT\{e,io) = (1 - e ) 



2n3/2 



(1- 



(20) 



(21) 



The chance of detecting the periastron precession in- 
creases with the geometrical factor GTv(e,^) which is re- 
lated to the alignment of the orbital ellipse with the line 
of sight. The optimal value for detecting the precession is 
lu = or 7r for small eccentricities, i.e. the semimajor axis 
should be perpendicular to the line of sight. For arbitrary ec- 
centricities, the optimal value for lu for the TTV observation 
is 



best 3 / 6e 

lutv ~ —it ± arccos . 

2 Vl + Vl + 24e 2 



(22) 
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which approaches 0° and 180° for small eccentricities as ex- 
pected, and 270° for large eccentricities. In case of e = 0.5, 
<j->extr = {235.4°, 304.6°}. The least favorable value for ui oc- 
curs when G T v(e, w) = 0, i.e. at w?^ rst = {90°, 270°},. 

2.3.2 Transit Duration Variations 

Now let us turn to the TDV effect. The observed duration of 
the mth transit can be calculated in the same way, namely 



Hn 



Ho + Fm, 



(23) 



where F is the shift in the transit duration per orbit. This 
factor is 



F = P H = 27rG DV (e, u, b)H-^-, 

J sec 



and 



Gov{e,uj, ft) 



1 — 2b e cos uj 
1 — ft 2 1 + e sin t 



(24) 



(25) 



The optimal orientation w D v' for detecting the TDV ef- 
fect for fixed e, ft, and Po can be found by maximizing 
HGov(e,u},b). The result is 



best 



-7r ± arccos 



4e 



1 + VI + 8e 2 



(26) 



and the worst orientation is at = {90°, 270°}, just 

like for the TTV case. Comparing w^v* an d ^dv* it is clear 
that the most favorable orientation in terms of the two ef- 
fects are similar, hence the chance of detecting the periastron 
motion through transit timing variations or transit duration 
variations is correlated. Both effects go away if the eccen- 
tricity is oriented parallel to the line of sight. We also note 
that for moderate values of e and small impact parameters, 
|Gtv| ~ |Gdv| which also implies that the most favorable 
geometry for detecting either TTVs or TDVs is similar. 



2.4 Error analysis 

Next we estimate the parameter measurement precision of 
the TTV and TDV effects for future observations. We con- 
sider the repeated observation of a particular transiting sys- 
tem over a total timespan T tot , measuring the transit timing 
T m and duration H m for each transit with respective errors 
cr(T) and a(H). (We discuss the specific values of o(T) and 
cr(H) for transit observations in Section |3.1[) . For simplic- 
ity, let us assume that these measurements are equidistant 
and in total N independent transits are observed, i.e. the 
mth transit is observed if m = 0, d, . . . , (N — l)d, where 
d = T tat /(NP ). 



2.4-1 Transit Timing Variations 



Using equation (|19[) . we can fit a second-order polynomial 
to these observations with unknown coefficients To, P obs and 
D by minimizing the merit function 

,2^2 



2 

Xtv — 



E 

m=0,d,...,(N-l)d 



T m - (T + P obs m + Dm 2 



a(T) 



(27) 



The minimization of the above function results in a linear set 
of equations in the parameters pi = {To, P bs, D}. Assuming 



Gaussian errors, the parameter estimation cova riance matr ix 
can be found from the Fisher matrix method (|Finnll 19921 ): 



{SptSpj} = (T )ij 

Here T is the Fisher matrix defined as 

i bt± a ar« d 



E 



m— O.d, . 
nfid 



.,(N-l)d 



2 (T) dpi d Pj 



(28) 



(29) 



where T m is the fiducial value of T m given by equation (|19[) . 
The marginalized expected squared parameter estimation 
error is given by the diagonal elements of the covariance 
error matrix cr' 2 (pi) = (T~ 1 )u. In particular, the resulting 
uncertainty of D becomes 

v / T80o-(r) 



a(D) 




4) 

2N 2 > Jn 



(30) 



Here, the first equality is valid for arbitrary N, while the sec- 
ond is its first order approximation for large N. The leadi ng 
order approximation is verified against IPress et al.l (|l992T l . 



2-4-2 Transit Duration Variations 

We can repeat the same calculations as above for the ob- 
servation of the half transit duration H m to measure the 
variation factor F. The merit function in this case 



2 

Xdv 



E 

=0,d,...,(iV- 



l)d 



H m - (Ho + Fm) 



a(H) 



(31) 



has to be minimized for the same set of observations. This 
minimization again leads to a linear set of equations in the 
parameters Ho and F. The Fisher matrix method in this 
case gives the uncertainty in F as 



a(F) 



12a (H) 



ds/N(N 2 - 1) 

VU T tot Y + 2N 2 J VTV • 



(32) 



Figure [2] shows the detection significance of the TDV 
measurement \F\/a(F), if the precession rate in F is given 
by the general relativistic formula, equation (|24[) . Here each 
transit is assumed to be measured (i.e. d — 1) with a 
precision o~(H) = 5 sec for a total observation time of 
Ttot = 4 years. These assumptions are realistic for the fu- 
ture Kepler mission (see § l3.1l below). Other parameters are 
M* = A/q, R t = Rq, ft — 0, and we averaged over the 
possible orientations of w. For other parameters, 



\F\/c(F) 
\Fq\/g(Fq) 



1 



d^T^¥ 

5 sec 



Rq 

' Ttot 
4yr 



3/2 



1/3 



(33) 



implying that the detection significance can be even better. 

Figure [2] clearly shows that the chances of detecting the 
precession effects through the TDV effect is encouraging. 
The detection significance of the general relativistic preces- 
sion of a transiting exoplanet with eccentricity e > 0.2 and 
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Figure 2. The significance of detecting the GR periastron preces- 
sion /'b /ao(l' ) through the TDV effect as a function of orbital 
period and eccentricity. The transit duration is assumed to be 
measured for 4 years each with a 5 sec error, for a Sun-like star 
on a non-inclined orbit. Increasing the mass or radius of the star, 
or the impact parameter increases the detection significance (see 
text). 



period P < 5 days, is typically over the 1-cr level. Gener- 
ally, equation (|30[) and equation (|32|l can be used directly 
to check what kind of observations are required to detect 
the precession of the periastron through the TTV or TDV 
methods, respectively. 



3 THE CASE OF XO-3b, HD 147506b, GJ436b 
AND HD 17156b 

As of this writing, four TEPs are know n with a non-zero 
eccen tricity within 3-cr, namely XO- 3 b lljohns-Krull et all 



l2007fl. HD 147506b jBakos et al 



Butler et al. 2004; Gillon et al 



l2007t) . GJ 436bb 
2007), and HD 17156bb 



Fischer et alj|2007l; iBarbieri et al.ll2007l ). The planet TrES 



1 (jAlonso et al. 2004) has also been reported as an object 
with non-zero eccentricity, however, it is zero within 2-cr thus 
we omit from our analysis. We note here that recently both 
GJ 436 and HD 17156 have been suggested to h ave another 
plane tary companions (see iRibas et al.l 120081 ; IShort et al.l 
l200ct ). The secular period of the periastron motion are de- 
termined by the mass of the star M*, the orbital period Po, 
and the eccentricity e, while the timing variation constant 
D is also affected by the actual argument of pericenter, uj. 
The transit duration variation factor F is affected indirectly 
by the geometrical ratio a/P* and directly by the impact 
parameter b. These parameters are summarized in the first 
seven columns of Table Q] for these four ETEP systems. The 
derived GR periastron precession period, P scc can be found 
in the 8th column of the table. 

In addition to the inevitable periastron precession 
caused by GR, there might be other sources of perturba- 
tions causing periastron precession. The last column gives 
the minimum mass to semimajor axis ratio of a hypothet- 
ical exterior perturber (e.g. a planet or an asteroid belt), 
in Earth mass units, which causes the same periastron pre- 
session rate as tha t caused by th e gene ral relativity. This 
estimate based on IPrice fc Rush! (1 19791 ), and is valid for 



Table 2. Transit timing variation factor (D, in seconds) and 
the number of transits (iV20y,3— <r,2scc) what should be detected 
almost uniformly in a 20 year long timespan, each with an error 
of 2 second to confirm the precession within 3-o\ 



System 


D (seconds) 




N20y, 3-rr, 2 sec 


HD 147506b 


-(5.9 ±0.7) 


io- 7 


6600 


XO-3b 


+ (4.3 ±0.6) 


io- 7 


1280 


GJ 436b 


+ (5.0 ±0.5) 


IO" 8 


44160 


HD 17156b 


-(6.9 + 0.8) 


10- 8 


97 ■ 10 6 



a,2 > 3a and for non-resonant cases. Note that the mini- 
mum mass of the perturber scales with the third power of 
the semimajor axis ratio to cause a comparable precession 
rate as GR. The numbers show that the precession caused 
by additional planets in the system, if present, can easily 
cause a larger precession rate than GR. In case of orbital 
resonances with exterior planets, the precession rate ca n be 
even larger (|Holman fc Murray! 120051 : lAgol et al.ll2005l ). To 
be conservative, we examine whether the precession rate can 
be measured to a precision better than that corresponding 
to GR. 

To obtain a high significance, k-a detection of the pe- 
riastron precession using transit timing variations, we need 
ka(D) ~ \D\. Using equation (|30p . the total number of tran- 
sits necessary to measure D with this precision is 



JVtv = 180 



a(T) 



k\D\ 



Ttot 



(34) 



Thus, the number of such required transit observations is 
extremely sensitive to the orbital period Po and the ob- 
servation timespan. Table [2] gives the corresponding val- 
ues of the transit timing variation factors for the known 
ETEP systems and this number of observations, assuming a 
Ttot = 20 year long observational timespan, timing precision 
of o(T) — 2 sec and 3-cr sensitivity of the GR periastron pre- 
cession level. The table shows the recently discovered XO-3b 
system is a promising candidate to detect the GR periastron 
precession through the TTV effect, while the other ETEP 
systems require unrealistically many observations. We note 
that if other perturbing planets are present in these systems 
and lead to a precession rate that is larger by a factor of 
10 than the GR precession rate, then the number of detec- 
tions (during the same T tot timespan) is lower by a factor 
of 100. It is interesting that the best candidate (by far) is 
XO-3b even though its eccentricity is not as large as that of 
HD 147506b or HD 17156b. 

Let us now turn to the observational constraints for the 
detection of transit duration variations. Using equation (|32[) . 
the total number of required observations within T to t is 



A+ 



12 



a(H) 



k\F\ 



Po 
Ttot 



(35) 



Note that Nov is not as sensitive to the Po/Ttot ratio as 
Atv, and implies that a smaller number of observations is 
typically necessary. 

In Table [3]we present the values of the transit duration 
variation factor, F, and the number of required observations 
to reach the same 3-cr confidence for detecting the GR peri- 
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Table 3. Transit duration variation factor (F, in seconds) and 
the number of transits (N2 y , 3— a, 2 sec) what should be detected 
almost uniformly in a 4 year long timespan, each with an error of 
2 sec to confirm the precession within 3-cx. 



Table 4. Uncertainties of the transit time and transit duration 
measurements for the four known ETEPs, assuming Sloan z'- 
band photometric data taken with a 1 sec cadence and 1 mmag 
photometric precision. 



System 


F (seconds) 




-^4y, 3-CT, 2scc 


System 


<ti,i(T) (sec) 


0-1,1 (H) (sec) 


HD 147506b 


-(1.8 ±0.3) 


10" 2 


20 


HD 147506b 


5.3 


4.8 


XO-3b 


-(5.4 ±0.6) 


10- 3 


70 


XO-3b 


6.9 


4.7 


GJ 436b 


-(4.1 ± 0.5) 


• KT 3 


85 


GJ 436b 


6.7 


8.4 


HD 17156b 


-(3.2 ± 0.2) 


■ 10~ 3 


8970 


HD 17156b 


5.4 


6.1 



astron precessiorjf). Here we assumed a shorter observation 
timespan, Ttot = 4 year (i.e. shorter compared to the 20 year 
long timespan necessary for the detection of the TTV ef- 
fect), the same timing precision of a(H) = 2 sec and the 
same level of detection, 3-er. The best known ETEP system 
for TDV detection is therefore HD 147506b, but the number 
of necessary observations is feasible for XO-3b and GJ 436b 
as well. 

3.1 Photometric detection 

The precision for measuring the transit timing and transit 
duration for a photometric observation can be estimated as 
follows. Since the time of the ingress (Ti) and the time of the 
egress (T E ) - i.e. when the center of the planet crosses the 
limb of the star inwards or outwards, respectively - defines 
both the time of the transit center and the half duration like 

T = i(T E + Ti), (36) 
H = ^(Te-Tj), (37) 

moreover Ti and Te can be treated as uncorrelated vari- 
ables, therefore the uncertainties of the transit time and half 
duration would be nearly the same, i.e. <r(T) ~ ^(H). We 
have estimated these uncertainties for the four distinct plan- 
ets using Monte-Carlo simulations by fitting transit light 
curves on mock data sets. We have used the observed plane- 
tary parameters as an input for these artificial light curves. 
The fit was performed as suming quadratic limb darkening 
(see iMandel fc Agol I2002T ) in the Sloan z band. The mock 
light curves were sampled with Ati = 1 sec cadence and an 
additional Gaussian noise of a\(vti) — lmmag was added. 
The resulting uncertainties, (Ti,i(T) and g\^{H) for the four 
planets are presented in Table [4] Since the depth of the four 
transits are nearly the same (see the appropriate normalized 
radii, p = R p /R+, all between 0.068 < p < 0.085), the un- 
certainties ai i i(T) and ai t i(H) are almost the same for the 
four cases. Using these normalized values, one can easily es- 
timate the uncertainties for arbitrary sampling cadence At 
and photometric precision a(m) using 

*V) ~ o lA T)^J^, (38) 
1 mmag V 1 sec 

3 Note that since the measurement of the TDV effect relies on 
fitting 2 parameters, instead of 3 parameters for the TTV effect, 
the 3-<t confidence corresponds to a higher confidence level for the 
TDV effect. 



a(H) « .u^^S (39) 
1 mmag V 1 sec 

For comparison, note that the expecte d photometric preci - 
sion of the Kepler space telescope (see lBorucki et aLll2007r ) 
is 1 mmag for observing a light curve of a bright, M v =8.8 
star with a 1 sec sampling cadence. Since the star XO-3 
has almost the same apparent magnitude, it is clear, the 
transit durations would be detected with an accuracy of 
<y(H) « 5 sec if this star was in the field of Kepler. There- 
fore, equation (|32[) and Table [3] shows that the transit du- 
ration variations would be detectable for HD 147506(b) or 
XO-3(b)-like systems due to the GR periastron precession, 
within 3-<t confidence with a Kepler-type mission within ap- 
proximately 3 or 4 years, respectively. 



4 SUMMARY 

The first four eccentric transiting exoplanetary systems have 
been discovered during 2007. The precession of an eccentric 
orbit causes variations both in the transit timings and tran- 
sit durations. We estimated the significance of measuring the 
corresponding observable effects compared to the inevitable 
precession rate of general relativity. We applied these calcu- 
lations to predict the significance of measuring the effect for 
the four known eccentric transiting planetary systems. Our 
calculations show that a space-borne telescope is adequate 
to detect the change in the transit durations to a high signif- 
icance better than the GR periastron precession rate within 
a 3 - 4 year timespan (in a continuous observing mode). The 
same kind of instruments would need more than a decade 
to detect the corresponding transit time variations to this 
sensitivity even for the most optimistic known system. 

The C 0R0T mission has already found two transiting 
planets (see lBarge et al.| [2008; Alo nso et al.ll200"8l ) and there 
are two known pla nets in the planned field-of-view of th e 
Kepler mission (see lO'Donovan et al.ll2006l : IPal et al.ll2008T l. 
Our results suggest that if an eccentric transiting planet is 
found in the Kepler or C0R0T field, these missions will be 
able to measure the periastron precession rate to a very high 
significance within their mission lifetime or with the support 
of ground-based observations on a longer time scale. This 
will provide an independent test of the theory of general 
relativity and will also be useful for testing for the presence 
of other planets in these systems. 
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APPENDIX A: MEAN LONGITUDE AT THE 
TRANSIT INSTANCES 

The derivation of equation {TJ goes as follows. According to 
Kepler's equation, E — e sin E = M — A — lo, one can write 



A = u> + E — esin_B. The only thing what is to be done is 
to calculate the eccentric anomaly E for the instance when 
the orbiting body intersect the semi-line with the argument 
angle ipo- The latter means that the true anomaly v of the 
body is v = ipo — w, by definition. The relation between the 
eccentric and true anomaly is 



E 

tan — = 
2 



- e v 

— tan-, 
■ e 2, 



which is equivalent with 

e + cos v 



cos E 



sinE 



1 + e cos v ' 
Vl — e 2 sin t 
1 + e cos v 



(Al) 

(A2) 
(A3) 



Using the addition theorem, the sine and cosine of the angle 
lo + E can be written as: 



cos(w + E) 



sin(u + E) 



e + cos(((9o — <j-0 
1 + e cos(<£>o — u) 



\/l — e 2 sin(((9o — <^) 
1 + e cos(<,9o — 
e + cos(i^o — lo) 



1 + e cos(<^o — ^) 



+ 



+ COS LO 



Vl — e 2 sin((^o — lo) 



(A4) 



(A5) 



1 + e cos(</?o — 
Thus, the mean longitude itself is going to be 
A = lo + E — es'mE — arg [cos(w + E),sm(LO + E)\ — 

(A6) 



1 + e cos v 

If both arguments of the above arg[-, •] function is multiplied 
by the always positive common denominator 1 + ecos(<^o — 
lo), one gets after some simplification: 



LO + E 



arg k + cos ipo + 
h + sin ipo 



h(k sin ifo — h cos tpo) 
k(k sin ipo — h cos tpo) 



1 + VT 



(A7) 



Putting all terms together and replacing the appropriate 
terms by e± — k sin ipo — hcos tpo, en = k cos tpo + hsin tpo and 
£ = 1 — Vl — e 2 , we get equation |l]) . The partial derivatives 
of equation ([l]) become 



dk 



h 



dX 

dh + 2 



2-1 
k 



-(l-£) 



h + (2 + en) sin Lpo 
k + (2 + eii) cos Lpo 



(A8) 
(A9) 



APPENDIX B: TANGENTIAL VELOCITY AND 
POSITION AT THE TRANSIT 

It is known from the theory of the two-body problem that 
the angular momentum of a body orbiting around a mass of 
GM = fi and having an orbit with the semimajor axis of a 
and eccentricity e is C = \J fia(l — e 2 ). Since C = r«tan for 
all points, the tangential velocity would be 



C 



f'ta 



V Ata(l - 



-1 + ecos(y>o — w) 
1 a(l-e 2 ) 



(Bl) 
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Using Kepler's Third Law, i.e. /i = n 2 a 3 , the above equation 
can be reordered to 



Wtan = an 



1 + e cos(<^o — 



For (/po = tt/2, equation (|B2|) becomes 
1 + e sin a; 



"tan = an- 



(B2) 



(B3) 
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